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1. Scalar and vector currents 

Form factors fo and /+ can be extracted from scalar and vector 3 -point correlators — see 
diagram in Fig. 1. The scalar current is a local, conserved current 

(K\s\D)=fr K (g 2 ) A ^^ (i.i) 

moc - m)s 

with S = ^F^F. Here q is the difference of the four momenta of the mesons, q = pr> — Pk, and the 
process D — > Klv is used as an example. One or both of the mesons are given a spatial momentum 
p using so called twisted boundary conditions (i.e. a phase at the boundary): 

&(x + ejL) = e 2 * e &(x) (1.2) 

where L is the size of the lattice. This gives the quark a momentum 2n0 /L. Note that 9 can be 
tuned to get a desired value for q 2 , e.g. q 2 = 0. We use different kinematical set-ups: kinematics 
A, where only one of the mesons has spatial momentum and the other one is at rest (in Fig. 1, s 
quark would have the momentum), and kinematics C, where both mesons have the same spatial 
momentum (in Fig. 1, the light quark would have the momentum). We have tested this method 
carefully, for example by checking that the speed of light is one (see Fig. 2), and that the amplitude 
of the meson correlator depends on the momentum like 1 / y/E. 
The vector current can be written as 

M 2 -M 2 



(K\V»\D)=f^(q*) 



+r(« 2 )^/. (i.3) 

9 



We have chosen to use = y^, a tasteless, spatial vector current. This has to be a 1-link current, 
if we have Goldstone mesons. We also tested a local, temporal vector current / with a non- 
Goldstone D s (YsTt) f° r A? -> Ul (denoted as V t in Fig. 6). We use MILC Nf = 2 + 1 lattices to 
do the calculations — see Table 1 for more details. 



scalar or vector current 

s f c 

K C Z> D 

< 1 ^ 

time T 



Figure 1: Diagram of the 3-point correlator setup. 



1.1 Fitting 

To extract the form factors we do a simultaneous least % 2 fit to both 3-point correlators and 
the corresponding 2-point meson correlators. For a given semileptonic decay, say D K,we also 
fit all q 2 values simultaneously. For the 2-point correlators the fit function is the usual sum of 
exponentials (with the oscillating states, as we are dealing with staggered quarks): e.g. for the D 
meson 

Q,(0 - J £( b f) 2 ( e ~ Eft + e ~ Ef{N, ~ t) ) -K^^-iyie-^ + e"^^-')). (1.4) 
j k 
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Figure 2: Test: Speed of light. 

The fit functions for the 3-point correlators have similar form: e.g. for D — » K we have 

C M (tJ) =IlA^e- £ ?^')-^^ ( e- £ H^)(-l)^' (1.5) 

j k j k 

-IL^V^e^^^-lJ' + ^^^V^e-^^-^-l) 7 , (1.6) 

j k j k 

where Aoo = b^b® (K\S\D) / (2y/M D E K ) gives the desired form factor fo at a given q 2 . We use three 
or four time separations T for the mesons (see Fig. 1), as the 3-point correlators are oscillating. 

2. Renormalization of the currents 

The scalar current is absolutely normalized (note the bare quark masses in Eq. (1.1) — the 
renormalization factors cancel). However, the vector current does need to be renormalized. We 
extract the renormalization factor Z from the symmetric vector current by demanding f+^ H (0) = 1 
for H = D, D s , ris, and T] c . The extracted Z factors agree over a range of momenta and different 
mesons for both charm-charm and charm-strange currents — see Fig. 3. This is essential, as we 
want to use Z to renormalize a charm- strange and a charm-light current. 

The local, temporal current V t is renormalized using fo at g^ ax extracted from the scalar cur- 
rent. Note that at g^ ax the temporal vector current gives the form factor fo directly, as the coeffi- 



ensemble 


size, L 3 x N t 


physical size 


# configs. 


# time sources 


mi 


coarse 


20 3 x 64 


w (2.4 fm) 3 


2259 


8 


ph m s /3.5 


fine 


28 3 x 96 


w (2.4 fm) 3 


1911 


4 


Pa m s /4.2 



Table 1: Details of the MILC 2+1 flavor lattice configurations used in this study: lattice size, number 
of configurations, number of time sources per configuration, and light valence quark mass m/ (compared 
to strange quark mass m s ). The mass values for HISQ valence light quarks are tuned to match the same 
goldstone pion mass as those for the asqtad sea light quarks. The HISQ valence s quark masses are tuned to 
the physical value — see [2] for more details. 
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Figure 3: Renormalization factor Z for the 1-link vector current. 
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Figure 4: Test of the symmetric scalar current. 



?max) 



cient multiplying / + vanishes (see Eq. 1.3). We can thus calculate /o(^max)' anc ^ set /(w(#n 

/0,,s(^max)- 

As yet another test we calculated the symmetric scalar current, (H\S\H ), for H = D s , r\ s , and 
rj c . At q 2 = 0, and when lattice spacing a — » 0, this is expected to give 



(H\S q \H)(q l = 0) 



2_n\- dm H 



(2.1) 



This is indeed the case, as can be seen in Fig. 4. 



3. Preliminary results: form factors fo and /_ 



+ 



Our preliminary results for the D —> n and D K semileptonic decay form factors fo and 
/+ are presented in Figs. 5, 6, and compared to experimental results in Section 3.2. We consider 
different mesons: For example, we calculate the charm-strange current 3-point correlator with 
different spectator quarks, light, strange and charm. This allows us to compare the form factors for 
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Figure 5: Preliminary results: form factors, charm to light decay. 
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Figure 6: Preliminary results: form factors, charm to strange decay. 



D —> K, D s — » ri s and r\ c —> D s . The semileptonic decay form factors are very insensitive to the 
spectator quark: D — » K and D s r\ s form factors are almost identical. However, if the spectator 
quark is as heavy as the charm quark, as in r\ c ^D s , the form factors do have a noticeably different 
shape. The insensitiveness of the form factors to the spectator quark is also seen in the light-charm 
current: D —> % and D s — >> K form factors have the same shape. One would expect to see similar 
behaviour in the corresponding B/B s form factors, i.e. basically no dependence on the spectator 
quark. 

3.1 z-expansion and continuum extrapolation 

It is convenient to transform the form factors to z-space to do the continuum extrapolation. 
This is done as follows: First we remove the poles from the form factors, 

- (i - 4r)f^ K ^ /rv) - (i - (3-D 
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Figure 7: Form factors in z-space. 



Here D to K is used as an example, so the poles are M^* and M^* . We then change from q to z, 




-q 2 + Vt^ 



= — , t + = (m D + m K ) 



(3.2) 



— note that we have taken to = in the standard transformation formula — and fit the lattice data 
as power series in z, 



We include terms up to z 4 . Note that bo = co, because fo(q 2 = 0) = f+(q 2 = 0). The lattice spacing 
dependence is very small, and the extrapolation to a = is shown in Fig. 7. 

3.2 Comparison with experimental results 

In Figs. 8 and 9 we compare our new results to earlier results as well as experimental results. 
Earlier results for D K form factor at q 2 = from HPQCD Collaboration are from [3], and the 
experimental results by CLEO Collaboration are from [4]. Both the value at q 2 = and the shape 
of the D —> K form factors fo and /+ agree very well with experimental results, where |V^| is 
calculated from the CKM matrix assuming unitarity. 

4. Conclusions 

We have presented here the first results from a very high precision study of D meson semilep- 
tonic decay vector form factors. We have looked at semileptonic decays of different mesons, and 
have shown that the form factors are very insensitive to the spectator quarks — we expect the same 
to be true for BIB S processes. Indeed our results indicate that B and B s semileptonic form factors 
at a given q 2 should differ by less than 5%, whereas their decay constants differ by approximately 
20%. QCD sum rules [5] predict SU(3) breaking effects at around 10% for both. We need to repeat 
our calculation with a lighter sea light quark mass to do a chiral extrapolation, and the simulations 



fr K (z) = E b n (a)z n , f^ K (z) = £ c n {a)z n . 



(3.3) 



n>0 n>0 



6 



The D to K and D to 71 semileptonic decay form factors from Lattice QCD 



Jonna Koponen 





HPQCD 2011 data h 


i 


i 




HPQCD 2011 h 


— )K — i 






HPQCD 2010 


— )K — i 






CLEO-C2009 h 


— )K — i 










* 










-> 


l < 















0.005 0.01 0.015 0.02 

a 2 (fm) 2 

Figure 8: Form factor at g 2 = 0, extrapolated to a — 0. 
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Figure 9: D ^ K form factors and experimental results from CLEO. 
are already underway. 
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